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ABSTRACT 

A  baroclinic  model  programmed  for  numerical  solution  is  used  to  in- 
vestigate the  effects  of  diabatic  heating  on  the  development  of  pressure 
waves  for  periods  up  to  24  hours.   Several  forms  of  the  model  are  studied 
with  results  being  obtained  for  various  wavelengths.   The  results  lead  to 
some  meteorological  improbabilities;  however,  the  relative  development  of 
the  pressure  wave  appears  valid  for  short  time  periods. 
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LIST  OF  SYMBOLS 
Variables  and  Constants  -  Roman  Alphabet 

A  -  D         amplitudes  of  various  sine  and  cosine  functions 

(-*.  Rossby  wave  speed,—  \J_  —  /& / \^ 

E  -  H         shorthand  equivalent  of  various  equations  as 

designated  in  text,  (appear  with  or  without  primes) 

f  the  Coriolis  parameter,  Qj}   .sinLp 

fQ  an  average  value  of  the  Coriolis  parameter 

^  acceleration  of  gravity 

i  unit  vector  in  X  direction 

(K  unit  vector  in  Z  direction 

k  wave  number,  2  77*  /L 

L  wave  length 

p  atmospheric  pressure 

Q  heating  term 

Rz  amplitude  of  the  pressure  wave 

t  time 

JJL  mean  zonal  flow  in  X  direction 

<x_  component  of  perturbation  velocity  in  X  direction 

%/  horizontal  velocity  vector 

\r  component  of  velocity  in  Y  direction 

x  coordinate  toward  east 

y  coordinate  toward  north 

z  height  of  isobaric  surface 
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Variables  and  Constants  -  Greek  Alphabet 

c<  coefficient  of  certain  terms,    2  +  46  ^  "  ^  J 

o  ' 

@  variation  of  the  Coriolis  parameter,   3f/d  y 

C  vertical  component  of  relative  vorticity,  =£~^7  $ 

£  phase  angle 

8  potential  temperature 

6  stability  factor,  =  -£-  -^  JaLS 

0  geopotential,  ^* 

^f  geographical  latitude 

t*J  total  derivative  of  pressure  with  respect  to  time 

/L  the  angular  speed  of  the  earth's  rotation 

"t~  positive  tolerance 

Subscripts  and  Superscripts 
L  vertical  index  denoting  a  pressure  level 

r\  index  denoting  the  nth  estimate  of  a  quantity 

Operators 
S7  two  dimensional  del  operator  on  a  constant  pressure 

surface 

_  a 

V  Laplacian  operator 

£\  finite  difference  operator  with  respect  to  p  and  t 

l(  )|  absolute  value  of  a  quantity 


1.    Introduction 

Most  numerical  prediction  methods  for  obtaining  prognostic  charts  of 
the  pressure  field  neglect  diabatic  heating  effects  under  the  assumption 
that  their  influence  is  negligible  for  short  period  prognosis. 

Under  certain  circumstances,  however,  the  exchange  of  sensible  heat 
between  an  air  mass  and  the  underlying  surface  may  be  considerable.   It  is 
therefore  the  purpose  of  this  paper  to  investigate  the  effects  of  such 
diabatic  exchange  on  the  development  of  pressure  waves  for  periods  up  to 
24  hours  using  a  numerical  model. 


2.   The  Model 

The  model  with  which  the  aforementioned  heating  effects  were  studied 
was  developed  from  the  geostrophic  relationship,  the  vorticity  equation, 
and  the  thermodynamic  equation, 

V=-f-XVjZS;  (1) 

il  +  V-Vff  +  5)=fo2u2  .  (2) 

J&    +  V'Vli.  +   5cj  =  -  _Q.  .  (3) 

apat  a?  t 

The  dependent  variables  are  expressed  as  the  sum  of  a  base  state  plus 
perturbation  quantity  (indicated  by  a  prime), 

0=  0CPji)  +  0'(X,^?^);  (4) 

V=  u  (?)<?  +  v'(*,^?/fc  );  (5) 

6J=  ^'CX^,?,*),  (6) 

Q=  QC?)-t  Q'(p)  um-KQCp)«0,  (7) 

S  =  S  +  S  '  (8) 

Equations  (2)  and  (3)  are  then  linearized  by  the  usual  method  of  assum- 
ing that  they  are  satisfied  by  the  basic  quantities  alone  and  by  the  total 
quantities  (sum  of  the  basic  plus  perturbation).   Subtracting  the  equations 
for  the  basic  flow  from  the  total  equations  and  neglecting  products  of 
perturbation  quantities  yields 

it  ^x         ax  a? 


Equation  (10)  is  further  modified  by  substituting  for  the  third 
term 

*££    =  -  £jl  oLLL  (id 

yielding, 

£j!  +U£jkl  -£4l£LU+  6U--JSL.  d2) 

dfdt  sxif     ax    otp  r 

A  diagnostic  equation  for  the  vertical  velocity  is  obtained  by 
differentiating  (9)  with  respect  to  p,  taking  the  Laplacian  of  (12),  and 
subtracting,  resulting  in 

Solutions  are  now  assumed  to  be  of  the  form 

Zvs^   =    A  fp,*)  sin  kx  +  3  <?,*)  £©*  ky  ,                              (14) 

u>*  =  C  Cfyt)  s\*  U%     +3>(p,  t)  co*  hx  ,  (15) 

Ql*    <?,  (Si*  kx+J,  )    +   QAfcw^X  +  ia  ).  (16) 

Next  equations  (14)  through  (16)  are  substituted  into  equations  (9) 
and  (13),  and  the  coefficients  of  ^Jniex  <ani£  cos  kx  equated  giving 

^A.  =   h  c  3-j^;  _^c.  ,  (i7> 

33.  =  -  K-c^a  -JLi^&,  (18) 


2  3Jl  £^B  -j±£  til  -J:  (%  s'* I-  Q, cos^ h   <19> 

fo   <*f>     ^    ^p    C 

^?v  "IT"  c  -^     ^   **   c 


Equations  (17)  through  (20)  are  then  solved  numerically  as  an  initial 
value  problem  as  outlined  by  Rosenthal  [lj  .    Given  initial  values  of 
A  and  B,  equations  (19)  and  (20)  are  solved  to  obtain  C  and  D.   These  in 
turn  are  substituted  into  (17)  and  (18)  to  obtain  prognostic  values  for 
A  and  B.   The  computational  cycle  is  continued  for  the  desired  time  in- 
terval. 


3.   Procedures 

Using  finite  differencing,  the  diagnostic  equations,  (19)  and  (20), 
can  be  approximated  at  the  odd  levels  as 

c«a""'ei  +  CM-,»KAp(if,  (UhrU)^-4akAPi(t  -3..,)  (19.) 

~ %±) -LCqj^^-Q,  cos $n)st  -h  ; 

to  r  c 


■3), 


+  ^4)j.fQtCMt  4-p,^ni;)«-  U*  .  (20a) 


<  a  3,5, 19 

At  the  even  levels  C  and  D  are  obtained  by  simple  averaging, 

c* r  Sbttfid     ^  a  Snl2^  *  (21) 

is   2,4, 20 

Boundary  conditions  are  that  C  and  D  at  the  top  and  bottom  levels 
are  zero,  implying  a  vanishing  vertical  velocity  there. 

The  prognostic  equations,  (If)  and  (JL8)  are  approximated  at  the  even 
levels  by 


(17a) 

(18a) 


■2&P 


At  the  odd  levels,  A  and  B  are  obtained  by  averaging, 
Al  =  Au>   +*-.,  ;  "3.  -  _Bu,+  3;-i  .  (22) 

»**    J  J  J  J  ♦  c  »  »  •  I  ..' 


*  A 


A  vertical  increment  of  50  mb  is  used,  affording  considerable  vertical 
resolution.   This  yields  a  21  point  vertical  grid  with  0  mb  as  level  1  and 
1000  mb  as  level  21,  as  shown  in  figure  1. 

Figure  1  -  Finite  Difference  Grid 
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The  values  of  C  and  D  are  obtained  by  a  simple  relaxation  scheme  as 
follows: 

(23) 


f.     =  _J 


) 


*--A 


5  =  «;+  ?U  ) 


(24) 
(25) 


^=3,5, 19,  and  £t     Ft  f   f  *  ~  O,   then 


(26) 


6-  /«?,  /*£ 3 

Equations  (17a)  and  (18a)  are  solved  by  an  iteration  technique  as 

follows, 
*/"*  (*  +  &*)**;  «)+££  CG.  Ul+G^C-t+^l)   J,  (28) 

(29) 


s^,)(-t+*t)=3L  W-tAi  C&/  (-6;+^*  rnj<f-t+A*)J 


£  a  2,4,  ...20,     h  «  1,2, 

Estimates  at  time  t  +  J^  t  are  indicated  by  the  superscript  n  while 
final  estimates  of  quantities  at  time  t  are  not  superscripted.   Initial 
estimates  for  G\      (t+At)  and  GL        (-t-f£tJ    are 

G-.(l)  (-t+A-k)  =   GL   C-O,  (30) 

G.;0Vt^t;  =  Cr,#  (4)a  (3i) 


The  iterative  process  is  continued,  recalculating  C  and  D  for  each 
new  value  of  the  superscript,  until 


A^'C-t-KA*)-  A'^\<*^t))     <L  <£,  (32) 

3fc-  "*'  C-fe-fA<;  -3/^V-fc^tJ  I  ^  <"■  (33) 


When  the  above  conditions  are  fulfilled,  the  new  prognostic  values 
of  A  and  B  are  then  used  to  calculate  corresponding  values  of  C  and  D. 
The  process  is  continued  for  the  desired  time  interval. 

The  amplitude  of  the  pressure  wave  and  phase  angle  can  then  be 
calculated  using  the  relationships,      * 

*'  =  nm  (  cos  kx  -  C  )  ,  (34) 


R  =  V  A  +  b\  <35) 

z  ' 

tan  <J>   »  A/B  ,  <36) 

2 

The  model  was  programmed  for  solution  by  the  CDC  1604  computer  #f  the 
U.  S.  Naval  Postgraduate  School.   The  time  step,  /\   t,  was  taken  as  30 
minutes,  and  the  tolerance,   <f"  ,  was  1  cm.   Solutions  for  the  amplitude 
and  phase  angle  were  obtained  for  every  6  hours  up  to  48  hours. 

The  heating  function  was  introduced  in  the  layers  below  500  mb  as  a 
function  of  pressure  in  a  manner  suggested  by  Haltiner  [2]  .   For  this 
study,    Qm    was  sat  equal  to  zero  and 

$*  "  1*  (VVoo)*  (37) 

"rVoo 

In  (37),   ^   represents  the  amount  of  heating  being  considered  in 
joules/gr  sac,  p  is  the  particular  pressure  level  under  consideration 
below  500  mb,  and  p  6-00  is  the  500-mb  pressure.   For  the  purpose  of  this 
investigation,  <£  ^  was  taken  as  4  x  10"^   joules/gr  sec  which  is  approxi- 
mately equivalent  to  heating  of  400  langlays/day.   Haltiner  and  Wang  f3j 
propose  that  the  exchange  of  heat  between  an  air  mass  and  underlying  surface 
ftygtyfefe  in  excess  of  800  langleys/daya   However,  it  is  felt  that  the  800f 
figure  represents  an  extreme,  and  that  the  value  of  400  langleys/day  is  more 
representative* 

The  model  is  such  that  the  heating  function  is  constant  in  time  and 
space  for  a  given  problem,  while  the  phase  angle  and  amplitude  of  the  pres- 
sure wave  are  set  initially  and  than  allowed  to  develop. 

Various  permutations  of  the  model  were  tested  with  different  phase 
angles  for  the  heating  function^  and  with  different  wave  lengths. 


In  all  cases,  &    was  set  equal  to  .5  x  10   cgs  units,  while  /3    was 

-13   -1    -1 
equal  to  1.6  x  10    cm   sec   .   The  wave  lengths  considered  were  2000, 

6000,  and  10000  kilometers. 


4.    Results. 

The  results  are  expressed  as  a  comparison  between  a  no  heat  case  and 
heating  introduced  with  various  phase  angles. 

a.   The  first  model  studied  had  linearized  vertical  shear  with 
Lia      =  40  m/sec  and  OC  =  Uf(\  -  Tj  \     ,  i   =  50,  100. ..1000.   The  initial 
value  of  B  at  all  levels  was  97.4  metersf  and  A  at  all  levels  was  zero. 
The  phase  angle  of  the  pressure  wave  was  therefore  zero  at  all  levels. 

With  the  initial  input  of  heat  in  phase  with  the  pressure  wave,  the 
development  at  the  end  of  a  24-hr  period  in  the  lower  layers  was  consider- 
able for  all  wavelengths,  (figure  2).    However,  in  the  layer  850-700  mb,  a 
decrease  in  amplitude  appears  for  the  very  long  wave,  10000  km.   The  long 
wave,  6000  km,  showed  a  similar  pattern  of  development  with  a  minimum  of 
development  in  the  same  layer,  850-700  mb. 

The  movement  of  the  pressure  wave,  as  deduced  from  the  phase  angle, 
shows  that  the  heating  effect  caused  an  increased  movement  of  the  long 
waves  in  the  lower  layers,  (figure  3).   Without  heating,  the  very  long 
wave  retrogressed  slowly  at  all  levels  throughout  the  period.   For  the 
initial  6-hour  period  with  heating,  it  retrogressed  slowly  at  all  levels, 
then  quite  rapidly  in  the  lower  levels.   As  seen  in  figure  3,  at  the  end 
of  the  24-hr  period,  the  retrogression  in  the  upper  levels  is  less  than 
without  heating,  while  in  the  lower  levels  it  is  considerably  more.   In 
24  hours,  the  long  wave  progressed  approximately  one-quarter  wavelength 
without  heating.   With  heating,  it  retrogressed  initially  in  the  upper 
layers,  while  progressing  in  the  lower.   After  14  hours  it  progressed  at 
all  levels,  so  that  after  24  hours,  the  net  effect  is  progression.   The 
short  wave  progressed  over  three-quarters  of  a  wavelength  without  heating 
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at  all  levels,  and  slightly  less  with  heating. 

Choosing  the  500-mb  level  as  representative  of  the  upper  levels,  a 
gradual  increase  in  amplitude  with  time  and  heating  was  noted  for  all 
wavelengths,  (figure  4). 

The  same  model  with  the  heating  function  initially  180  degrees  out 
of  phase  with  the  pressure  wave,  showed  that  at  the  end  of  24  hours,  the 
two  waves  increased  in  amplitude  below  750  mb,  while  all  wavelengths 
increased  below  900  mb,  (figure  5).   The  10000-km  wave  was  the  only  one, 
when  compared  to  the  no  heat  case,  which  showed  a  larger  amplitude  at  all 
levels.   The  6000-km  wave  exhibited  a  marked  decrease  at  all  upper  levels, 
while  the  short  wave  decreased  slightly. 

The  very  long  wave  retrogressed  in  all  but  the  lowest  layers.   The 
short  wave  increased  its  progressive  movement  to  approximately  one  wave- 
length, while  the  long  wave  decreased  its  movement,  progressing  slowly  in 
the  lower  layers,  and  retrogressing  in  the  upper  layers,  (figure  6). 

A  look  at  figure  7  for  the  development  at  500-mb  clarifies  the  case 
of  the  very  long  wave.   It  too  decreased  in  amplitude  until  after  approxi- 
mately 18  hours  had  elapsed.  This  was  true  for  all  other  levels. 

For  initial  phase  angles  of  the  heating  function  of  90  and  270  degrees, 
the  same  marked  increase  in  amplitude  in  the  lower  levels  is  evident,  (fig- 
ures 8  and  9).   In  the  upper  levels,  for  an  initial  angle  of  90  degrees, 
the  short  wave  decreased  nearly  uniformly,  the  long  waves  increased.   For 
an  initial  angle  of  270  degrees,  the  opposite  appeared. 

b.   A  second  variation  of  the  basic  model  included  a  linear 
vertical  shear  to  500  mb,  and  then  a  parabolic  variation  of  the  zonal 
wind  to  the  top  of  the  atmosphere.   In  this  case,  &-soo      equaled  20  m/sec, 
and  &-3SS   was  30  m/sec.   A  further  modification  was  to  allow  for  tilt  of 
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the  pressure  wave  by  varying  the  phase  angle  in  accordance  with  the 
relationship,  phase  angle  equals  ~7T/3,   ( — t-  )    L  s   q    <o /ooo. 

"PlOOO   i 

In  all  cases,  this  model  showed  extreme  development  below  900  mb  for 
all  wavelengths,  regardless  of  the  initial  heating  angle. 

When  the  heating  phase  angle  was  zero,  the  10000-km  wave  developed 
very  slightly  at  upper  levels,  and  decreased  in  the  layer  650-800  mb,  (fig- 
ure 10).   The  6000-km  wave  developed  well  above  700  mb,  but  with  a  minimum 
of  development  at  800  mb.   The  short  wave  evinced  a  nearly  uniform  increase 
in  development  at  all  levels  above  800  mb. 

For  heating  with  an  initial  180-degree  phase  angle,  the  2000-km  wave 
showed  slight  additional  development  above  700  mb  except  for  a  slight  de- 
crease in  the  layer  300-200  mb,  and  the  10000-km  case  showed  increased 
development  throughout,  (figure  11).   Again,  an  examination  of  the  500- 
mb  level  reveals  that  the  very  long  wave  did  decrease  in  amplitude  until 
beyond  12  hours,  (figure  12). 

A  270-degree  heating  phase  angle  produced  increased  amplification  for 
all  levels  for  the  short  wave,  while  both  long  waves  decreased  above  800 
mb,  (figure  13). 

A  90-degree  heating  phase  angle  resulted  in  an  increased  amplitude  of 
the  6000-km  wave  at  all  levels,  (figure  14).   The  10000-km  wave  showed  a 
decreased  amplitude  from  900=800  mb,  and  increased  above  800  mb.   The 
short  wave  had  decreased  development  at  all  levels  above  900  mb. 

c.   A  number  of  other  cases  were  studied,  however,  the  results 
were  comparable  to  those  reported,  and  therefore  are  not  shown. 
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5.   Conclusions. 

The  restrictions  imposed  upon  the  basic  model  lead  to  certain 
meteorological  improbabilities,  i.e.,  excessive  speed  of  movement  in 
certain  cases  and  over  amplification  after  twenty-four  hours.   However, 
it  is  felt  that  for  short  time  periods,  the  relative  development  in  the 
layers  above  900  tab  is  plausible  even  on  a  quantitative  basis. 
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APPENDIX  I 
Figures  2  through  14  referred  to  in  text  are  included  here  for 
easy  reference. 
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